REM  ---- Exziteration --------------------------------
rem   
REM      f(z)=c*exp(z),    c=a+i*b
rem
rem       初期値を与えての反復計算
rem       初期値を与えないときは，z=0 が初期値
REM       nl<= n <=nu の範囲の数列をプリント
REM      |f(z)|>10^25 or Re(Xn)>60  で発散とみなす。
REM------------------------------- 5.18 '17, Joe -------
rem
PRINT"初期値を与えて Zn+1=c*exp(Zn) の反復を計算"
INPUT PROMPT"C=a+i*b をinputせよ，a,b=":a,b
line input prompt"初期値をあたえるか？　y/n?":fa$
if fa$="N" or fa$="n" then goto 10
input prompt"x0,y0=":x0,y0
goto 20
10 LET  x0=0
   LET  y0=0
20 LET  r=1.E+25
   print"表示する数列の最小と最大の n をinputせよ"
   input prompt"nl,nu=":nl,nu
   rem
   rem     計算
   rem
   for n=1 to nu
      LET  x1=EXP(x0)*(a*COS(y0)-b*SIN(y0))
      LET  y1=EXP(x0)*(b*COS(y0)+a*SIN(y0))
      LET  d=SQR(x1*x1+y1*y1)
      if nl<=n and n<=nu then goto 30 else goto 35
30    print n;")";x1;" , ";y1
35    IF x1>60 OR d>r THEN GOTO 40 ELSE GOTO 45
40    print"発散です d=";d;" n=";n
      goto 100
45    LET  x0=x1
      LET  y0=y1
   next n
   rem
100 END